Search for: All records

Abstract Globular clusters (GCs) provide statistically significant coeval populations of stars spanning various evolutionary stages, allowing robust constraints on stellar evolution model parameters and ages. We analyze eight old Milky Way GCs with metallicities between [Fe/H] = −2.31 and −0.77 by comparing theoretical isochrone sets from the Dartmouth Stellar Evolution Program to Hubble Space Telescope (HST) observations. The theoretical isochrones include uncertainties introduced by 21 stellar evolution parameters such as convective mixing, opacity, diffusion, and nuclear reactions, capturing much of the quantifiable physics used in our code. For each isochrone, we construct simulated color–magnitude diagrams (CMDs) near the main-sequence turnoff region and apply two full-CMD-fitting methods to fit HST Advanced Camera for Surveys data across a range of distances and reddening and measure the absolute age of each GC from the resulting posterior distribution, which accounts for uncertainties in the stellar models, observations, and fitting method. The resulting best-fitting absolute ages range from ≈11.5 to 13.5 Gyr, with a typical error of 0.5–0.75 Gyr; the data show a clear trend toward older ages at lower metallicities. Notably, distance and reddening account for over 50% of the uncertainty in age determination in each case, with metallicity,αabundance, mixing length, and helium diffusion being the most important stellar physics parameters for the error budget. We also provide an absolute age–metallicity relation for Milky Way GCs. 

Abstract We investigate the impact of massive primordial black holes (PBHs;m_BH ∼ 10⁶M_⊙) on the star formation and first galaxy assembly process using high-resolution hydrodynamical simulations fromz= 1100 toz ∼ 9. We find that PBH accretion is self-regulated by feedback, suppressing mass growth unless feedback is weak. PBHs accelerate structure formation by seeding dark matter (DM) halos and gravitationally attracting gas, but strong feedback can delay cooling and suppress star formation. In addition, the presence of baryon-DM streaming creates an offset between the PBH location and the peaks induced in gas density, promoting earlier and more efficient star formation compared to standard ΛCDM. Byz ∼ 10, PBH-seeded galaxies form dense star clusters, with PBH-to-stellar mass ratios comparable to observed high-zactive galactic nuclei like UHZ-1. Our results support PBHs as viable supermassive black hole (SMBH) seeds but do not exclude alternative scenarios. We emphasize that PBH-seeding provides a natural explanation for some of the newly discovered overmassive SMBHs at high redshift, in particular those with extreme ratios of BH-to-dynamical (virial) mass that challenge standard formation channels. Future studies with ultra-deep JWST surveys, the Roman Space Telescope, and radio surveys with facilities such as the Square Kilometre Array and Hydrogen Epoch of Reionization Array will be critical in distinguishing PBH-driven SMBH growth from other pathways. 

ABSTRACT In the local Universe, star formation is typically inefficient both globally and when considered as the fraction of gas converted into stars per local free-fall time. An important exception to this inefficiency is regions of high gravitational accelerations g, or equivalently surface densities $$\Sigma = g/(\pi \, G)$$, where stellar feedback is insufficient to overcome the self-gravity of dense gas clouds. In this paper, I explore whether dark matter can play an analogous role in providing the requisite accelerations on the scale of entire galaxies in the early cosmos. The key insight is that characteristic accelerations in dark matter haloes scale as $(1+z)^2$ at fixed halo mass. I show this is sufficient to make dark matter the source of intense accelerations that might induce efficient star formation on galactic scales at cosmic dawn in sufficiently massive haloes. The mass characterizing this regime scales as $$(1+z)^{-6}$$ and corresponds to a relatively constant comoving number density of $$n(>\!M_{\rm {vir}}) \approx 10^{-4}\, {\rm Mpc}^{-3}$$ at $$z \gtrsim 8$$. For somewhat rarer haloes, this model predicts stellar masses of $$M_{\star }\sim 10^{9}\, {\rm M}_{\odot }$$ can form in regions that end up with sizes $$\mathcal {O}(100\, {\rm pc})$$ over $$40\, {\rm Myr}$$ time-scales at $$z\approx 12-14$$; these numbers compare well to measurements for some of the brightest galaxies at that epoch from JWST observations. Dark matter and standard cosmological evolution may therefore be crucial for explaining the surprisingly high levels of star formation in the early Universe revealed by JWST. 

ABSTRACT JWST has revealed a large population of UV-bright galaxies at $$z\gtrsim 10$$ and possibly overly massive galaxies at $$z\gtrsim 7$$, challenging standard galaxy formation models in the ΛCDM cosmology. We use an empirical galaxy formation model to explore the potential of alleviating these tensions through an Early Dark Energy (EDE) model, originally proposed to solve the Hubble tension. Our benchmark model demonstrates excellent agreement with the UV luminosity functions (UVLFs) at $$4\lesssim z \lesssim 10$$ in both ΛCDM and EDE cosmologies. In the EDE cosmology, the UVLF measurements at $$z\simeq 12$$ based on spectroscopically confirmed galaxies (eight galaxies at $$z\simeq 11\!-\!13.5$$) exhibit no tension with the benchmark model. Photometric constraints at $$12 \lesssim z\lesssim 16$$ can be fully explained within EDE via either moderately increased star-formation efficiencies ($$\epsilon _{\ast}\sim 3\!-\!10\ \hbox{per cent}$$ at $$M_{\rm halo}\sim 10^{10.5}{\, \rm M_\odot }$$) or enhanced UV variabilities ($$\sigma _{\rm UV}\sim 0.8\!-\!1.3$$ mag at $$M_{\rm halo}\sim 10^{10.5}{\, \rm M_\odot }$$) that are within the scatter of hydrodynamical simulation predictions. A similar agreement is difficult to achieve in $$\Lambda$$CDM, especially at $$z\gtrsim 14$$, where the required $$\sigma _{\rm UV}$$ exceeds the maximum value seen in simulations. Furthermore, the implausibly large cosmic stellar mass densities inferred from some JWST observations are no longer in tension with cosmology when the EDE is considered. Our findings highlight EDE as an intriguing unified solution to a fundamental problem in cosmology and the recent tensions raised by JWST observations. Data at the highest redshifts reached by JWST will be crucial for differentiating modified galaxy formation physics from new cosmological physics. 

We use FIRE-2 zoom simulations of Milky Way size disk galaxies to derive easy-to-use relationships between the observed circular speed of the Galaxy at the Solar location,v_c, and dark matter properties of relevance for direct detection experiments: the dark matter density, the dark matter velocity dispersion, and the speed distribution of dark matter particles near the Solar location. We find that both the local dark matter density and 3D velocity dispersion follow tight power laws withv_c. Using this relation together with the observed circular speed of the Milky Way at the Solar radius, we infer the local dark matter density and velocity dispersion near the Sun to beρ= 0.42±0.06 GeV cm^-3and^σ_3D= 280⁺¹⁹_-18km s^-1. We also find that the distribution of dark matter particle speeds is well-described by a modified Maxwellian with two shape parameters, both of which correlate with the observedv_c. We use that modified Maxwellian to predict the speed distribution of dark matter near the Sun and find that it peaks at a most probable speed of 257 km s^-1and begins to truncate sharply above 470 km s^-1. This peak speed is somewhat higher than expected from the standard halo model, and the truncation occurs well below the formal escape speed to infinity, with fewer very-high-speed particles than assumed in the standard halo model. 

ABSTRACT Fuzzy dark matter (FDM), comprised of ultralight ($$m \sim 10^{-22}\,{\rm eV}$$) boson particles, has received significant attention as a viable alternative to cold dark matter (CDM), as it approximates CDM on large scales ($${\gtrsim}1$$ Mpc) while potentially resolving some of its small-scale problems via kiloparsec-scale quantum interference. However, the most basic FDM model, with one free parameter (the boson mass), is subject to a tension: small boson masses yield the desired cores of dwarf galaxies but underpredict structure in the Lyman-α forest, while large boson masses render FDM effectively identical to CDM. This Catch-22 problem may be alleviated by considering an axion-like particle with attractive particle self-interactions. We simulate an idealized FDM halo with self-interactions parametrized by an energy decay constant $$f \sim 10^{15}~\rm {GeV}$$ related to the axion symmetry-breaking conjectured to solve the strong-CP problem in particle physics. We observe solitons, a hallmark of FDM, condensing within a broader halo envelope, and find that the density profile and soliton mass depend on self-interaction strength. We propose generalized formulae to extend those from previous works to include self-interactions. We also investigate a critical mass threshold predicted for strong interactions at which the soliton collapses into a compact, unresolved state. We find that the collapse happens quickly, and its effects are initially contained to the central region of the halo. 

Abstract We analyze the first cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modeled as a constant fractional energy dissipation (f_diss= 0.75). In this paper, we examine the properties of dwarf galaxies withM_*∼ 10⁵–10⁹M_⊙in both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm²g⁻¹. On the contrary, models with (σ/m) = 10 cm²g⁻¹produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm²g⁻¹are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the Hirotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with (σ/m) ≳ 0.1 cm²g⁻¹,f_diss= 0.75 are in tension in massive dwarfs (M_halo∼ 10¹¹M_⊙) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures. 

ABSTRACT Recent observations with JWST have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the big bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox$$^{\it HR}$$, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments (FIRE) project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox$$^{\it HR}$$ re-simulates the cosmic volume ($L=22.1$ cMpc) of the original FIREbox run with eight times higher mass resolution ($$m_{\rm b}\sim {}7800\, M_\odot$$), but with identical physics, down to $$z\sim {}6$$. FIREbox$$^{\it HR}$$ predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at $$z\sim {}6{\!-\!}14$$, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox$$^{\it HR}$$, the SFE–halo mass relation for intermediate mass haloes ($$M_{\rm halo}\sim {}10^9{\!-\!}10^{11}\, {\rm M}_\odot$$) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE–halo mass relation lead to a larger contribution from lower mass haloes at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE–halo mass relation inferred from FIREbox$$^{\it HR}$$ allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at $$z\gt 12$$ will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn. 

Abstract We present ∼300 stellar metallicity measurements in two faint M31 dwarf galaxies, Andromeda XVI (M_V= −7.5) and Andromeda XXVIII (M_V= –8.8), derived using metallicity-sensitive calcium H and K narrowband Hubble Space Telescope imaging. These are the first individual stellar metallicities in And XVI (95 stars). Our And XXVIII sample (191 stars) is a factor of ∼15 increase over literature metallicities. For And XVI, we measure $〈\mathrm{[Fe/H]}〉=-{2.17}_{-0.05}^{+0.05}$, ${\sigma }_{\mathrm{[Fe/H]}}={0.33}_{-0.07}^{+0.07}$, and ∇_[Fe/H]= −0.23 ± 0.15 dex $R_e^{-1}$. We find that And XVI is more metal-rich than Milky Way ultrafaint dwarf galaxies of similar luminosity, which may be a result of its unusually extended star formation history. For And XXVIII, we measure $〈\mathrm{[Fe/H]}〉=-{1.95}_{-0.04}^{+0.04}$, ${\sigma }_{\mathrm{[Fe/H]}}={0.34}_{-0.05}^{+0.05}$, and ∇_[Fe/H]= −0.46 ± 0.10 dex $R_e^{-1}$, placing it on the dwarf galaxy mass–metallicity relation. Neither galaxy has a metallicity distribution function (MDF) with an abrupt metal-rich truncation, suggesting that star formation fell off gradually. The stellar metallicity gradient measurements are among the first for faint (L≲ 10⁶L_⊙) galaxies outside the Milky Way halo. Both galaxies’ gradients are consistent with predictions from the FIRE simulations, where an age–gradient strength relationship is the observational consequence of stellar feedback that produces dark matter cores. We include a catalog for community spectroscopic follow-up, including 19 extremely metal-poor ([Fe/H] < –3.0) star candidates, which make up 7% of And XVI’s MDF and 6% of And XXVIII’s. 

ABSTRACT We introduce an analytic surface density profile for dark matter haloes that accurately reproduces the structure of simulated haloes of mass Mvir = 107–1011 M⊙, making it useful for modelling line-of-sight (LOS) perturbers in strong gravitational lensing models. The two-parameter function has an analytic deflection potential and is more accurate than the projected Navarro, Frenk, and White profile commonly adopted at this mass scale for perturbers, especially at the small radii of most relevant for lensing perturbations. Using a characteristic radius, R−1, where the log slope of surface density is equal to −1, and an associated surface density, Σ−1, we can represent the expected lensing signal from LOS haloes statistically, for an ensemble of halo orientations, using a distribution of projected concentration parameters, $$\mathcal {C}_{\rm vir} := r_{\rm vir}/ R_{-1}$$. Though an individual halo can have a projected concentration that varies with orientation with respect to the observer, the range of projected concentrations correlates with the usual three-dimensional halo concentration in a way that enables ease of use.